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Abstract 

Using a lattice-based Monte Carlo code for simulating self-avoiding flexible polymers in three 
dimensions in the absence of explicit hydrodynamics, we study their Rouse modes. For self-avoiding 
polymers, the Rouse modes are not expected to be statistically independent; nevertheless, we 
demonstrate that numerically these modes maintain a high degree of statistical independence. Based 
on high-precision simulation data we put forward an approximate analytical expression for the mode 
amplitude correlation functions for long polymers. From this, we derive analytically and confirm 
numerically several scaling properties for self-avoiding flexible polymers, such as (i) the real-space 
end-to-end distance, (ii) the end-to-end vector correlation function, (iii) the correlation function of 
the small spatial vector connecting two nearby monomers at the middle of a polymer, and (iv) the 
anomalous dynamics of the middle monomer. Importantly, expanding on our recent work on the 
theory of polymer translocation, we also demonstrate that the anomalous dynamics of the middle 
monomer can be obtained from the forces it experiences, by the use of the fluctuation-dissipation 
theorem. 

PACS numbers: 36.20.-r,64.70.km,82.35.Lr 
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I. INTRODUCTION 



Polymer dynamics is a field where first-principle analytical derivations of most quantities 
of interest from microscopic monomeric movements are difficult to come by. Indeed, the 
answer to the question why phantom polymers (theoretical realizations of polymers that 
can intersect themselves) remain, to this day, the central pillar for studying polymer dy- 
namics is easily traced to the fact that they allow full analytical calculations of essentially 
all of their dynamical properties. In the absence of explicit hydrodynamics, a phantom 
polymer's dynamics is described by the so-called Rouse equation that forms the basis of 
all analytical calculations pj, 12|. The Rouse equation holds in the high- viscosity limit of 
the surrounding medium, and the dynamics of the monomers are "overdamped", i.e., the 
velocity of each monomer at any given time is proportional to the total force it experiences. 
In the Rouse equation the total force in any monomer is comprised of the spring forces 
due to its neighboring monomers, and the random thermal forces from the surrounding 
medium. The linearity of the Rouse equation allows one to decompose it into linear dynam- 
ical equations of independently evolving (Rouse) modes, which are the Fourier transforms 
of the monomer co-ordinates in three dimensional space. In detail, we consider a polymer 
of length A^, consisting of + 1 monomers connected sequentially by N bonds (harmonic 
springs). The positions -R„(t) of monomers n = 0, . . . , at time t can then undergo Fourier 
transformation, yielding for the amplitude of the p-th mode {p is an integer > 0) 

N 



n=0 



7r(n + l/2)p 

ivTi 



Rn{t). (1) 



The cornerstone of all analytical calculations for phantom polymer dynamics is the following 
relation for p,q ^ 0, derived exactly from the Rouse equation: 



- - N 
Xp,(t) = {Xp{t) . X,(0)) oc - exp 



im< (2) 



where 5pq is the Kronecker delta function, and A a constant. All throughout this paper, the 
angular brackets represent an average over the equilibrium ensemble of polymers. Equation 
(El) is further supplemented by X^pit) = for p 7^ 0, and Xoo(t) = m6{t), where Xo(t) is 
the location of the center-of-mass of the polymer at time t, and m oc 1/A^ is the mobility 
of the center-of-mass of the polymer. The Kronecker delta terms signify the statistical 
independence of the modes of a phantom polymer. Using these mode amplitude correlation 



functions, the quantities of interest for a phantom polymer can be analytically tracked by 
reconstructing them from the modes |2|. However, in reality, a polymer is not phantom, 
but is self-avoiding. The self-avoidance property introduces long-range correlations along 
the backbone of the polymer, and it destroys the linearity of the Rouse equation, making 
a similar [to Eq. ([2])] equation for the mode amplitude correlation functions impossible to 
derive from the appropriately formulated Rouse equation. Because of this reason, it comes as 
no surprise to us that for self-avoiding polymers we have not been able to find, in published 
literature, a comprehensive study of the mode amplitude correlations that connects to the 
scaling properties of self-avoiding polymers. 

The purpose of this paper is to put forward an approximate analytical expression for 
the mode amplitude correlation functions for long polymers, which is then corroborated 
with extensive simulations of a lattice-based model for the dynamics of self-avoiding fiexible 
polymers in three dimensions in the absence of explicit hydrodynamics. To be more precise, 
with u ^ 0.588 (the Flory exponent in three dimensions), and Ai and A2 two constants, we 
demonstrate that for a self-avoiding Rouse polymer of length N 



X,,it) = (Xp(t) ■ X,(0)) ^ A, exp 



5pg (3) 



for p, g 7^ holds up to a very good approximation for long polymers, while Xop(t) = for 
p 7^ and Xoo(t) oc {1/N)6(t) holds exactly. The reasons why Eq. ([3]) is not exact are the 
following: (a) although numerically the modes maintain a high degree of statistical indepen- 
dence, the mode amplitude correlation functions are not exactly statistically independent 
(we do not expect to be so anyway); and (b) for higher modes {p, g > 6) there are small devi- 
ations from exponential behavior at long times. Nevertheless, given that Xpp{0) is a rapidly 
decaying function of p, the dominant contribution of the modes to quantities of interest for 
the polymer at long time-scales come from the lower modes. As a result, such quantities can 
be analytically reconstructed from Eq. ([3]). We demonstrate this by using Eq. (I3|) to derive 
several scaling properties for self-avoiding polymers, such as (i) the real-space end-to-end 
distance, (ii) the end-to-end vector correlation function, (iii) the correlation function of the 
small spatial vector connecting two nearby monomers at the middle of the polymer, and 
(iv) the anomalous dynamics of the middle monomer. Some of these scaling laws can also 
be obtained from the dynamic scaling law p|. Note that except the case for characterizing 
the anomalous dynamics of the middle monomer, all the other quantities (as above) concern 



only the polymer's internal structure; implying that for (i-iii) we only need Eq. ((31), while for 
the anomalous dynamics of the middle monomer we also need Xoo(t) and Xop(t) for p 7^ 0. 
Note also that the standard Rouse result Eq. ([2]) is recovered when we simply replace v by 
1/2 for a phantom polymer. Apparently, the dominant consequence of volume exclusion is 
the different size scaling of the self-avoiding chains. The inability for polymers to cross each 
other does not seem to cause large cross-correlations between different modes. Thus, given 
the content of this paper, we expect that Eq. (I3|) will provide researchers a way forward for 
analytical treatment of the properties of self-avoiding polymers. 

The structure of this paper is as follows. In Sec. [Til we describe our polymer model and 
demonstrate the scaling property Eq. (jSj). In Sec. Uni we use Eq. ([3]) to obtain the scal- 
ing properties for the real-space end-to-end distance and the end-to-end vector correlation 
function. In Sec. [TVl we derive the scaling properties of the small spatial vector connect- 
ing two nearby monomers. In Sec. [Vl we derive the anomalous dynamics of the middle 
monomer, and by expanding our recent work on the theory of polymer translocation, show 
that the anomalous dynamics of the middle monomer can also be obtained by using the 
fluctuation-dissipation theorem. Our conclusions are then summarized in Sec. |VT1 



II. OUR POLYMER MODEL AND THE SCALING PROPERTIES OF THE MODE 
AMPLITUDE CORRELATION FUNCTIONS 

A. Our polymer model and the calculation of the mode amplitudes 

Over the last years, we have developed a highly efficient simulation approach to polymer 
dynamics. This is made possible via a lattice polymer model, based on Rubinstein's repton 
model for a single reptating polymer, with the addition of sideways moves (Rouse dy- 
namics). A detailed description of this model, its computationally efficient irnplementation 
and a study of some of its properties and applications can be found in Refs. jj, Ol. 

In this model, each polymer is represented by a sequential string of monomers, living on a 
face-centered-cubic lattice with periodic boundary conditions in all three spatial directions. 
Monomers adjacent in the string are located either in the same, or in neighboring lattice 
sites. The polymers are self-avoiding: multiple occupation of lattice sites is not allowed, 
except for a set of adjacent monomers. The polymers move through a sequence of random 
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single-monomer hops to neighboring lattice sites. These hops can be along the contour 
of the polymer, thus explicitly providing reptation dynamics. They can also change the 
contour "sideways", providing Rouse dynamics. Each kind of movement is attempted with 
a statistical rate of unity, which provides us with the definition of time. This model has 
been used before to simulate the diffusion and exchange of polymers in an equilibrated layer 
of adsorbed polymers Q|. Recently, we have used this code extensively to study polymer 



translocation under a variety of circumstances 



dynamics of polymer adsorption 
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12| . and also to study the 



Given that our model has periodic boundary conditions in all three spatial directions, we 
use the following definition for the mode amplitude: 

'7r{n + l/2)p 



N + 



1 ^ 

n=0 



Ro{t)] cos 



N+1 



(4) 



where Ro{t) is the location of the zero-th monomer at time t. In this definition we only 
need the spanning vectors, i.e., monomer co-ordinates w.r.t. that of monomer zero. These 
spanning vectors are obtained from a summation over the bond vectors between monomers 
adjacent in the string. This completely avoids the invocation of periodic images, and allows 
for spanning distances which exceed half the box size. Note that Eq. (Jlj) can be derived 
from Eq. ([T]) as follows. First we use 



N + 



Ecos 

n=0 



n{n + l/2)p 



cos 



n{n + l/2)p' 
N + 1 



1a 

2 "p>p'' 



(5) 



in order to express Rn{t), the location of the n-th monomer w.r.t. that of the center-of-mass 
Remit) = Xo(t), from the inverse Fourier transform of Eq. JT]), as 

'TT{n + l/2)p 



Then we write 



Rn{t) = Remit) + 2 cos 
p=l 



Roit) = Remit) + 2^Xp(t) cos 
p=l 



N+1 



7rp 
2(iV+l) 



I.e., 



i?,(t)-i?o(t) = 2j]Xp(t) <jcos 
p=i 



nin + l/2)p 



cos 



Tip 



(6) 



(7) 



(8) 



A^ + 1 J "■"'[2(iV + l 

from which Eq. (Jlj) follows with the use of Eq. (I5|) . Equation ((41) shows that the modes g 7^ 
relate only to the polymers' structural configuration, as remarked earlier. 



B. The scaling properties of the mode amphtude correlation functions 



We start with the behavior of Xoo(t). We note that the internal forces over the entire 
polymer at any time sum to zero, and therefore, (for the overdamped dynamics in the 
Rouse model) the motion of its center-of-mass is simply proportional to the average thermal 
fluctuating force on the entire polymer. This force is 5-correlated in time, which implies that 
the center-of-mass of the polymer simply performs a random walk. In fact, Xoo(t) = m'6{t), 
and the mobility m' oc This leads us to the result that Xoo(t) oc {1/N)6{t). We will 

return to the motion of the center-of-mass in Sec. |Vl 
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FIG. 1: (a) The scaling Xpp(O) ~ N^'^/p^+^'' demonstrated for p = 1, . . . , 10 for iV = 200, 400, 600 
and 1000. The solid line corresponds to NXpp{0) ~ {N/p)-^^+'^''^ w (iV/p)-2-i76_ (t.) xhe matrix 



Xpq{0) = Xpq{0) / ^y Xpp{0)Cqq{0) in logarithmic grayscale for p,q = 1, . . . , 10 for = 1000 in a 
square checkerboard plot. See text for details. 

Next, for the behavior of Xop(t) for p 7^ we proceed as follows. We note that Xop(O) = 0. 
Indeed, Xop(O) is the equilibrium average over the dot product of Xq and Xp at zero time 
difference, for which, for any given value of Xp, the center-of-mass is equally likely to be 
present at any location within the periodic box. For the calculation of Xop(t), having 
summed over all the locations of Xq for a given value of Xp at the first step, and then having 
summed over the configurations Xp at the second step, we obtain Xop(O) = (from the first 
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step). Thereafter, we use the result that the center-of-mass performs a simple random walk, 
further implying that Xop(t) = Xop(O), which equals zero. 




(p/N) t 



FIG. 2: Data collapse for -\n[Xpp{t) / Xpp{{))] as a function of {p/NY+'^''t, for p = 1 (stars), 2 
(triangles down), 3 (crosses), 4 (triangles up), 5 (diamonds), 6 (pluses), 7 (circles), and N = 200 
(red), 400 (blue), 600 (green) and 1000 (orange). The collapsed data are fitted well with an 
exponential; however, we note that there are small (but systematic) deviations from the exponential 
behavior for p » 1. The sohd (black) line corresponds to -\n[Xpp{t) / Xpp{<d)] ~ [p/XY^'^'^t. 



Finally, we demonstrate the scaling property of Xpq{t) for g 7^ in Figs. [T] and [2l 
In Fig. [11(a) we present the data for p = 1, . . . , 10 for = 200, 400, 600 and 1000, 
and demonstrate that Xpp(O) ~ N"^^ 1^^+"^^ up to a high precision. In Fig. [I](b) we plot 
the matrix Xpg(O) = -^pg(0)/A/^CT'(0)^?g(0) logarithmic grayscale for p, g = 1, . . . , 10 
for = 1000, wherein the diagonal elements are unity by construction. The off-diagonal 
elements of Xpq{^) ^tre typically three or more orders of magnitude smaller than the diagonal 
ones. Cross-correlations between two even or between two odd modes are strictly not zero. 
Cross-correlations between even and odd modes are much smaller, and it is likely that these 
are also strictly not zero, but we cannot ascertain that within our numerical precision. We 
emphasize that the off-diagonal elements of Xpg(O) not being zero is not caused by the lack 
of numerical precision, or is due to artifacts of our model — similar features have been found 
in Ref. {l^ for single polymers, and in Ref. 1^ for polymer melts; the modes are simply not 



statistically independent for self-avoiding polymers. Nevertheless, the precision with which 
the modes remain statistically independent is remarkable. 

In Fig. M we plot - \n[Xpp{t) / Xpp{0)] as a function of {p/Ny+^^t for p = 1, . . . , 7, and 
N = 200, 400, 600 and 1000, and obtain a data collapse. The collapsed data are well fitted 
by an exponential behavior of Xpp(t)/Xpp{0)] however, we note that there are small (but 
systematic) deviations from the exponential behavior for p > 6. These deviations are larger 
for larger p, although this fact is not very clearly discernible in Fig. [2l 

As remarked before Eq. ([3]), Figs. [T] and [2] collectively demonstrate that the scaling be- 
havior Eq. ([3]) is not exact; instead, it is a rather good approximation. In the following 
sections we will use the approximate scaling behavior of Eq. ([3]) and demonstrate that it 
reproduces the scaling behavior of several observables for self-avoiding polymers. 



III. SCALING OF THE END-TO-END DISTANCE AND THE END-TO-END 
VECTOR CORRELATION FUNCTION 



In order to calculate the equilibrium end-to-end distance we use Eq. ([8]) and write 



RM{t)-Ro{t) = 2j2Xp{t) <jcos 



7r{N + l/2)p 



N + 1 

from which, using the scaling relation ([3j) we extract 

p e odd 



— cos 



Tip 



2(iV+l) 



(9) 



irNp 
_2{N + l) 



p G odd 



P^ 



sm 



np 



2{N+1)_ 



(10) 



On the r.h.s. of Eq. (fTOl) the first term p sums up to a numerical constant, and 

p G odd 

the second term [7rp/2(A^ + 1)] can be converted to an integral 

p G odd 



8Ai / dx 



sin^[7ra;/2] 

'0 

with an integrable singularity at x = 0. The first term produces the well-known Flory scaling 
behavior {[Riy(t) — _Ro(^)]^) ~ N"^" , while the second term provides a correction of 0(1) to 
the scaling which can be neglected in the scaling limit. 
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The end-to-end vector correlation function can be similarly expressed in terms of the 
mode amplitude correlation functions; however, in order to obtain its scaling behavior we 
do not need to do so. The point is that at long times the dominant contribution of the end- 
to-end vector comes from the mode p = 1, and thus the approximate scaling ([3]) confirms 
that the end-to-end vector correlation function decays exponentially in time, for which the 
characteristic time for the decay is given by the Rouse time tr ~ A^^^^'^. 

IV. THE CORRELATION FUNCTION OF THE SMALL SPATIAL VECTOR 
CONNECTING TWO NEARBY MONOMERS AT THE MIDDLE OF THE POLY- 
MER 



Consider the small spatial vector r„(t) from monomer {N — n)/2 to monomer {N + n)/2 
with n N. In this section we study the behavior of C„(t) = (^n(^) ■ ^n(O)). 
By expanding in terms of modes, we have 



p=i 



-4 Mt) 

p g odd 



sm 



7r{N + n + l)p' 
2{N+1) 

71 np 



— cos 



7r(A^ - n + l)p 
2{N + 1) 



_2{N + l] 

Thereafter, using the approximate scaling ([3]) we obtain 



(11) 



sm 



8Ain 



2u 



p £ odd 

°° dx 



imp 



2{N + 1] 



exp 



-A, 



p 



l+2u 



2u 







l+2u 



X 



sin'[7rx/2] exp[-A2x'^^''it/n'-^^''^ 



(12) 



We did not find a closed-form expression for the integral in Eq. ( fT2l ): hence, in order to 
obtain the scaling behavior of (r„(t) ■ r„(0)) with time, we proceed with an approximation 
in the following manner. First, we note that at long times the dominant contribution to the 
integral comes from x < A2{t/n^'^'^'^y^^^^'^''\ This observation allows us to write 



C„(t) ^ 8A,n 



2u 



^2 ^l/(l+2iy) 



dx 

„l+2u 



sin^[7ra;/2]. 



(13) 



'0 ^ 

Next, for such (small) values of x in the integral (fT3l ). we can approximate sin(7ra;/2) by 
7rx/2, leading to 

42^1/(1+2,.) 



C„(t) ^ 2Ai'K^n 



X 



l-2v 



2n'A^Af-''^n'^"^-'h~'^. 



(14) 



The behavior C„(t) ~ t-2(i-!-)/(i+2i/) ^^^^ ^^^^^ ^j^g Rouse time, the lifetime of the 
lowest mode {p = 1) in the summation of Eq. (fT2ll . In Fig. [3] we numerically confirm this 
behavior of C„(t) for N = 390 and n = 6. 
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FIG. 3: Plot of Cn{t) for N = 390 and n = 6. The solid line corresponds to C„(t) ~ 
^-2{i-u)/(i+2u) ^-0.378 rjij^g steep drop in the data at very long times corresponds to the terminal 
(exponential) Rouse relaxation. 



V. ANOMALOUS DYNAMICS OF THE MIDDLE MONOMER 



A. Mean-square displacement of the middle monomer 



To obtain the mean-square displacement of the middle monomer we first use Eq. ^ to 
write 

= RN/2{t) - RN/2iO) = Xoo(t) - Xoo(O) + 2 J2 [^p(^) - ^p(O)]- (15) 

p G even 

We then use Xop(t) = and the approximate scaling relation ([3]) to obtain 

Xpp{t) 

p G even '- 



^ {[Remit) - RcmiO)?) + § Xpp(O) 



Xpp(O) 



([i?,„(t) - Rcmm^) + ^AiN^" |l - exp 



p S even 



l+2u 1 



-A 



_P 



(16) 
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Given that the center-of-mass performs a random walk, the first (center-of-mass) term on 
the r.h.s. of Eq. (fT6l) . as argued in Sec. [Tll increases hnearly with t, while in the limit 
N oo the second term can be converted to an integral: 



p 6 even 



-A. 



P 



l+2v 



2u 



8Ai 



dx 



[1 - exp(-A2a;^+2'^t)] 



dx 



/ -TT^ [1 - eM-A2x'^ni (17) 

'0 



which, up to some constant factor, scales with time as ti+a-^. Once again, the ^1+2" behavior 
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FIG. 4: The scaling behaviors for the mean-square displacement (Ar^(t)): for the center-of-mass 
(bottom set of points), and for the middle monomer (top set of points). The dot-dashed line 
corresponds to (Ar^(f)) ~ t and the solid line corresponds to (Ar^(t)) ~ iiu/{i+2u) ^ 

of the mean-square displacement of the middle monomer will only hold till the lifetime of 
the lowest mode in the summation of Eq. (fTTll {p = 2) scaling as the Rouse time tji ~ A^^"^^*^, 
after which {Ar'^^^) to increase linearly with t. Note that Eqs. (fTBl) and (fTTll confirm 
the well-known result that by the Rouse time, the middle monomer typically displaces itself 
by the spatial extent of the polymer ~ N^" . 

In Fig. m we confirm the above behavior of the mean-square displacement of the middle 
monomer and that of the center-of-mass. We have also checked, in support of Xoo(t) ~ S{t), 
that the probability distribution of the displacement of the center-of-mass of the polymer 
along all three spatial direction is Gaussian, with the width scaling ~ \/t (data not shown). 

11 



B. Fluctuation dissipation theorem and the anomalous dynamics of the middle 
monomer 



In Sees. [TTl and [V] we argued that the force on the center-of-mass of the polymer is the 
average over all the thermal forces on the monomers, leading us to Xoo(t) ~ S(t), which 
yielded the ~ t behavior of the mean-square displacement of the center-of-mass. This raises 
the natural question, namely: is it possible to derive the anomalous dynamics of the middle 
monomer — i.e., the (Ar^^g) ~ up to the Rouse time and ~ t thereafter — from the 
combined forces (internal and external) that act on the middle monomer? 

We trace back the formulation of this problem to our recent work on unbiased poly- 

n 

mer translocation through a narrow pore in a membrane [8|. Therein we showed, in the 
following manner, that the anomalous dynamics of translocation — essentially that of the 
translocating monomer — can indeed be derived from the forces it experiences. The velocity 
of translocation v{t) (along the direction perpendicular to the membrane) is related to the 
force 0(t) (also acting upon it along the direction perpendicular to the membrane), via the 
"impedance" memory kernel fi{t) and "admittance" memory kernel a{t), by 

0(t) = f dt' n{t-t')v{t') + g{t), (18) 
Jo 

and 

v{t)= [ dt' a{t-t')(p{t') + h{t). (19) 
Jo 

In Eqs. (fTSl) and (fTOl) . g{t) and h{t) are the noise terms satisfying {g{t)) = {h(t)) = 0, and 
the corresponding fluctuation-dissipation theorems {g{t)g{t')) = (0(t)0(t'))t,=o = \f^(t — t')\ 
and {h(t)h{t')) = {v(t)v{t'))^=o = \a(t — t')\. Moreover, the uniqueness of the relation 
between v(t) and (f){t) dictates that their Laplace transforms must satisfy the condition 
fi{z)a{z) = 1. Finally, we characterized the anomalous dynamics of translocation by inte- 
grating (f (t)f (t'))0=o twice in time: we obtained the mean-square displacement (s^(t)) of 
the translocating monomer to behave as until the Rouse time tr, beyond which (s^(t)) 
linearly with t. We showed that the memory kernel scheme works beautifully for translo- 



cation out of planar confinements [9||, translocation by a pulling force at the head of the 
polymer (and additional back-pulling voltage jlli), field-driven translocation 12 1, and 



polymer adsorption jl3l |. 
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It is easy to see in the above analysis that if (0(t)0(t'))t,=o = — t')\ ~ for 
some a, then the mean-square displacement of the translocating monomer has to increase 
as t". We now show that the above scheme established for the anomalous dynamics of 
polymer translocation also works, with the modification that we now have to deal with 
vector velocity v{t) of the middle monomer and vector force f{t) on it, for the anomalous 
dynamics of the middle monomer of a self-avoiding polymer; namely that we will show that 
\fi{t — t')\ ~ t~i+2^^ which would imply (Ar^^g) ~ . Of course such power-law behavior 
would only hold till the Rouse time. 

In order to do so, we first argue, following Eq. f fTSll . that the impedance memory kernel 
fi(t) — that connects the forces on and the velocities of the middle monomer — scales 
~ t 1+2"' . Consider the thought experiment in which we grab the middle monomer and 
move it by a small distance 6r and hold it at its new position; this corresponds to v(t) = 
SrS{t). In time, the "information" that the middle monomer has moved to a new position 

at t = will propagate along the backbone of the polymer, and at time t, all the monomers 

1 

within a backbone distance rit ~ tw<^ — following the Rouse scaling — will equilibrate to 
this new situation. These rit equilibrated monomers are however stretched by an amount 
Sr. With the entropic spring constant of n equilibrated monomers scaling ~ n~'^'^ , the 
(restoring) force the middle monomer would experience at its new position is given by 
fit) ~ "(—Sr) ~ t~i+2^(— (5r) [force = (spring constant) x (stretching distance)]. In 
other words, /i(t) ~ t i+Z'' . The fiuctuation-dissipation theorem then dictates that (/(t) • 

/(t'))k-=o = W-t')\ ~ (t-t')"^ and (t/(t)-tT(t'))l/=o = ~ (t-t')^^- Note 

once again that these relations only hold till the Rouse time; by the Rouse time the entire 
polymer is equilibrated to the new position of the middle monomer, and beyond that time 
the forces on the middle monomer are uncorrelated. Finally, the anomalous dynamics of the 
middle monomer — (Ar^j^^) ~ t^^ up to the Rouse time and ~ t thereafter — is retrieved 
by integrating {v(t) ■ 'y(t'))/=o twice in time. 

Implementing the above thought experiment into practice and thereby tracking the restor- 
ing force on the middle monomer in order to determine fi(t) is a complicated task. Instead, 
we focus on the relation {f [t) ■ f [t'))^=o ~ t~^^ . We fix the position of the middle monomer 
of a self-avoiding polymer (this corresponds to v = 0) and take snapshots of it at equal in- 
tervals of time. We then take the snapshot at time t, evolve the entire polymer — with its 
middle monomer free — over At = 1 unit of time for a multiple number of times, and obtain 
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the average displacement vector u(t) of the middle monomer. Note that At = 1 is small 
enough such that the corresponding Eq. (fTOl) — with 0(t) and v{t) replaced by f{t) and 
v{t) respectively, and {h(t)) = — shows that u(t) oc f(t)y=o. The behavior of ('u(t) ■ uit')) 
in time then yields the behavior of (/(t) ■ /(t'))k=o = \f^{t ^ Determined in the above 
manner, we confirm \^(t)\ ~ t~i+2^ exp(— t/r^) in Fig. [5l with 2z//(l + 2z/) ^ 0.54, from 
which, as described in the above paragraph, the scaling (Ar^^g) ~ can be derived. 

10"' I — — — I 




t 



FIG. 5: The scaling behavior of for N = 400, determined in the manner described in the 

text. The sohd (black) curve corresponds to = 0.275t~°-^'' exp(-t/Tij), with tr 223801 for 

N = 400, obtained from best-fitting the corresponding data of Fig. [2l 



VI. CONCLUSION 

To conclude, in this paper we put forward an approximate analytical expression ([3]) for 
the mode amplitude correlation functions Xpq{t) for long polymers, and corroborate its 
accuracy using a lattice-based Monte Carlo simulations of self-avoiding flexible polymers in 
three dimensions in the absence of explicit hydrodynamics. We report that (a) the mode 
amplitude correlation functions are not exactly statistically independent (we do not expect 
them to be so anyway); instead, numerically the modes maintain a high degree of statistical 
independence, and (b) for higher modes {p, g > 6) there are small deviations from the 
exponential behavior ([3]) at long times. Nevertheless, given that Xpp{0) is a rapidly decaying 
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function of p as per Eq. ([3]), the dominant contribution of the modes to quantities of interest 
for the polymer at long time-scales come from the lower modes. As a result, such quantities 
can be analytically reconstructed from Eq. ((Sj). We demonstrate this by using Eq. (JS]) to 
derive several scaling properties for self-avoiding polymers, such as (i) the real-space end- 
to-end distance, (ii) the end-to-end vector correlation function, (iii) the correlation function 
of the small spatial vector connecting two nearby monomers at the middle of the polymer, 
and (iv) the anomalous dynamics of the middle monomer. Given the content of this paper, 
we expect that Eq. ([3]) will provide researchers a way forward for analytical treatment of 
the properties of self-avoiding polymers. 

Importantly, expanding on our recent work on the theory of polymer translocation, we 
also demonstrate that the anomalous dynamics of the middle monomer can be obtained 
from the forces it experiences, by the use of the fluctuation-dissipation theorem. Our hope 
is that in cases where a polymer's anomalous dynamics exponents are difficult to identify, 
use of the fluctuation-dissipation theorem will prove to be an indispensably useful tool; flrst 
such cases were our recent work on the theory of polymer translocation as well as above. 

In experimental situations, the polymer dynamics is often dominated by hydrodynamic 
interactions. In future work, we therefore intend to extend this study to polymer models 
with explicit hydrodynamics. 
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